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Abstract. 

Cosmic rays have been observed up to energies 10* times larger than those of the 
best particle accelerators. Studies of astrophysical particles (hadrons, neutrinos and 
photons) at their highest observed energies have implications for fundamental physics 
as well as astrophysics. Thus, the cosmic high energy frontier is the nexiis to new 
particle physics. This overview discusses recent advances being made in the physics 
and astrophysics of cosmic rays and cosmic 7-rays at the highest observed energies 
as well as the related physics and astrophysics of very high energy cosmic neutrinos. 
These topics touch on questions of grand unification, violations of Lorentz invariance 
as well as Planck scale physics and quantum gravity. 



1. Introduction 

Owing to the uncertainty principle, it has long been realized that the higher the particle 
energy attained, the smaller the scale of physics which can be probed. Thus, optical, 
UV and X-ray observations led to the understanding of the structure of the atom, 7-ray 
observations led to an understanding of the structure of the atomic nucleus and deep 
inelastic scattering experiments with high energy electrons led to an understanding of 
the structure of the proton. Accelerator experiments have led to an understanding 
of QCD and it is hoped that they will eventually reveal new physics which has been 
predicted at the TeV scale. To go beyond this, to questions of grand unification, and 
even Planck scale physics, one must turn to the cosmic accelerators with which Nature 
has provided us. As we will see, manifestations of physics with large extra dimensions 
can also be searched for at ultrahigh energies. 

Observations have been made of ultrahigh energy cosmic rays up to 300 EeV 
(3 X 10^ TeV) and of cosmic 7-rays up to 50 TeV energy. As of this writing, no very 
high or ultrahigh energy cosmic neutrinos have been detected, however, the AMANDA 
(Antarctic Muon and Neutrino Detector Array) experiment, now in operation, is 
searching for neutrinos above 1 TeV energy (Wischnewski 2002). We will review here the 
present status of the observations of various cosmic particles at their highest observed 
energies and the implications of these observations for astrophysics as well as physics. 
We will take a synoptic view of ultrahigh energy hadrons, photons and neutrinos. In 
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Figure 1. The ultrahigh energy cosmic ray spectral data from the analysis of Fly's 
Eye (closed triangles), AGASA (closed circles), HiRes I-monocular (open circles), and 
HiRes Il-monocular (open triangles) observations. 



this way, one can gain insights into the profound connections between different fields of 
observational astronomy and astrophysics which use different experimental techniques. 

2. The Highest Energy Cosmic Rays 

2.1. The Data 

Figure 1 shows the data (as of this writing) on the ultrahigh energy cosmic ray spectrum 
from the Fly's Eye, AGASA and HiRes detectors. Other data from Haverah Park and 
Yakutsk may be found in the review by Nagano and Watson (2000) are consistent with 
Figure 1. The new HiRes data are from Abu-Zayyad et al. (2002). 

For air showers produced by primaries of energies in the 1 to 3 EeV range, 
Hayashida, et al. (1999) have found a marked directional anisotropy with a A.5a excess 
from the galactic center region, a 3.9a excess from the Cygnus region of the galaxy, 
and a 4. Oct deficit from the galactic anticenter region. This is strong evidence that EeV 
cosmic rays are of galactic origin. A smaller galactic plane enhancement in EeV events 
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Figure 2. Average depth of shower maximum (Xmax) vs. energy compared to the 
calculated values for protons (upper curves) and Fe primaries (lower curves) (from 
Gaisser 2000; see references therein). 



was also reported by the Fly's Eye group (Dai et al. 1999). 

As shown in Figure 2 , at EeV energies, the primary particles appear to have a 
mixed or heavy composition, trending toward a light composition in the higher energy 
range around 30 EeV (Bird, et al. 1993; Abu-Zayyad, et al. 2000). This trend, together 
with evidence of a flattening in the cosmic ray spectrum in the 3 to 10 EeV energy range 
(Bird, et al. 1994; Takeda et al. 1998) is evidence for a new component of cosmic rays 
dominating above 10 EeV energy. 

The apparent isotropy (no galactic-plane enhancement) of cosmic rays above 10 
EeV {e.g., Takeda, et al. 1999), together with the difficulty of confining protons in the 
galaxy at 10 to 30 EeV energies, provide significant reasons to believe that the cosmic- 
ray component above 10 EeV is extragalactic in origin. As can be seen from Figure 1, 
this extragalactic component appears to extend to an energy of 300 EeV. Extention of 
this spectrum to higher energies is conceivable because such cosmic rays, if they exist, 
would be too rare to have been seen with present detectors. We will see in the next 
section that the existence of 300 EeV cosmic rays gives us a new mystery to solve. 
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Figure 3. The GZK cutoff energy, defined as the minimum energy predicted for a 
flux decrease of 1/e owing to intergalactic photomeson production interactions, as a 
function of redshift (Scully and Stecker 2002). 



2.2. The GZK Effect 

Thirty eight years ago, Penzias and Wilson (1965) reported the discovery of the cosmic 
2.7K thermal blackbody radiation which was produced very early on in the history of 
the universe and which led to the undisputed acceptance of the "big bang" theory of the 
origin of the universe. Much more recently, the COBE (Cosmic Background Explorer) 
satellite confirmed this discovery, showing that the cosmic background radiation (CBR) 
has the spectrum of the most perfect thermal blackbody known to man. COBE data 
also showed that this radiation (on angular scales > 7°) was isotropic to a part in 10^ 
(Mather et al. 1994). The perfect thermal character and smoothness of the CBR proved 
conclusively that this radiation is indeed cosmological and that, at the present time, it 
fills the entire universe with a 2.725 K thermal spectrum of radio to far-infrared photons 
with a density of ~ 400 cm~^. 

Shortly after the discovery of the CBR, Greisen (1966) and Zatscpin and Kuz'min 
(1966) predicted that pion-producing interactions of ultrahigh energy cosmic ray protons 
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with CBR photons of target density ~ 400 cm~^ should produce a cutoff in their 
spectrum at energies greater than ~ 50 EeV. This predicted effect has since become 
known as the GZK (Greisen-Zatsepin-Kuz'min) effect. Following the GZK papers, 
Stacker (1968) utilized data on the energy dependence of the photomeson production 
cross sections and inelasticities to calculate the mean energy loss time for protons 
propagating through the CBR in intergalactic space as a function of energy. Based 
on his results, Stecker (1968) then suggested that the particles of energy above the GZK 
cutoff energy (hereafter referred to as trans-GZK particles) must be coming from within 
the "Local Supercluster" of which we are a part and which is centered on the Virgo 
Cluster of galaxies. Thus, the "GZK cutoff" is not a true cutoff, but a suppression of 
the ultrahigh energy cosmic ray flux owing to a limitation of the propagation distance 
to a few tens of Mpc. 

The actual position of the GZK cutoff can differ from the 50 EeV predicted by 
Greisen (1966). In fact, there could actually be an enhancement at or near this energy 
owing to a "pileup" of cosmic rays starting out at higher energies and crowding up 
in energy space at or below the predicted cutoff energy (Puget, Stecker and Bredkamp 
1976; Hill and Schramm 1985; Berezinsky and Grigor'eva 1988; Stecker 1989; Stecker and 
Salamon 1999). The existence and intensity of this predicted pileup depends critially on 
the flatness and extent of the source spectrum, {i.e., the number of cosmic rays starting 
out at higher energies), but if its existence is conflrmed in the future by more sensitive 
detectors, it would be evidence for the GZK effect. 

Scully and Stecker (2002) have determined the GZK energy, defined as the energy 
for a fiux decrease of 1/e, as a function of redshift. At high redshifts, the target photon 
density increases by (l+z)^ and both the photon and initial cosmic ray energies increase 
by {1 + z). The results obtained by Scully and Stecker are shown in Figure 3. 

Figure 4 gives further results of Scully and Stecker (2002) compared with the present 
data. It shows the form of the cosmic ray spectrum to be expected from sources with 
a uniform redshift distribution and sources which follow the star formation rate. The 
required normalization and spectral index determine the energy requirements of any 
cosmological sources which are invoked to explain the observations. Pileup effects and 
GZK cutoffs are evident in the theoretical curves in this figure. As can be seen in 
Figure 4, the present data appear to be statistically consistent with either the presence 
or absence of a pileup effect. Future data with much better statistics are required to 
determine such a spectral structure. 

The AG AS A collaboration has recently reevaluated their energy determination and 
they claim a significant number of events at trans-GZK energies (Takeda, et al. 2003). 
However the analysis of observations made with the HiRes monocular detector array 
show only one event signiflcantly above 100 EeV and appear to be consistent with the 
GZK effect (Abu-Zayyad et al. 2002; see Figure 4 and section 3.7.) In fact, De Marco et 
al. (2003) have argued that the discrepency between the AG AS A and HiRes results is 
not statistically significant and that many more events are needed in order to determine 
whether or not there is a GZK effect. 



1e+26 



1e+25 ■ 




1e+23 ■ 



1e+22 



1e+17 



1e+18 



1e+19 
Energy (eV) 



1e+20 



1e+21 



Figure 4. Predicted spectra for cosmic ray protons as compared with the data. The 
middle curve and lowest curve assume an £;^2.75 gguj-ce spectrum with a uniform 
source distribution and one that follows the z distribution of the star formation rate 
respectively. The upper curve is for an £;~2.35 source spectrum which requires an order 
of magnitude more energy input and exhibits a "pileup effect" . 

The significance of a non-observation of a GZK effect is profound. Such a result 
either requires a large overdensity of UHECRs within about 100 Mpc, which would 
have to be emitted by "local" sources and trapped by magnetic fields, or it requires new 
physics such as the violation of Lorcntz invariance at ultrahigh energies. We will discuss 
these points in further detail below. 

2.3. Acceleration and Zevatrons: The "Bottom Up" Scenario 

The apparent lack of a GZK cutoff (with the exception of the new HiRes results) has 
led theorists to go on a hunt for nearby "zevatrons", i.e., astrophysical sources which 
can accelerate particles to energies 0(1 ZeV = lO^^eV). 

In most theoretical work in cosmic ray astrophysics, it is generally assumed that 
the diffusive shock acceleration process is the most likely mechanism for accelerating 
particles to high energy. (See, e.g., Jones (2000) and references therein.) In this case. 
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Figure 5. A "Hillas Plot" showing potential astrophysical zevatrons (from Olinto 
2000). The lines are for B vs. L for Emax = 0.1 ZeV for protons and iron nuclei as 
indicated. 

the maximum obtainable energy is given by Emax = keZ{u/c)BL, where u < c is the 
shock speed, eZ is the charge of the particle being accelerated, B is the magnetic field 
strength, L is the size of the accelerating region and the numerical parameter k — 0{1) 
(Drury 1994). Taking k — 1 and u — c, one finds 

E^ax = 0.9Z{BL) 

with E in EeV, B in /iG and L in kpc. This assumes that particles can be accelerated 
efficiently up until the moment when they can no longer be contained by the source, 
i.e. until their gyroradius becomes larger than the size of the source. Hillas (1984) used 
this relation to construct a plot of B vs. L for various candidate astrophysical objects. 
A "Hillas plot" of this kind, recently constructed by Olinto (2000), is shown in Figure 
5. 

Given the relationship between Emax and BL as shown in Figure 5., there are 
not too many astrophysical candidates for zevatrons. Of these, galactic sources such 
as white dwarfs, neutron stars, pulsars, and magnetars can be ruled out because their 
galactic distribution would lead to anisotropics above 10 EeV which would be similar 

to those observed at lower energies by Hayashida et al. (1999), and this is not the 
case. Perhaps the most promising potential zevatrons are radio lobes of strong radio 
galaxies (Biermann and Strittmattcr 1987). The trick is that such sources need to be 
found close enough to avoid the GZK cutoff {e.g., Elbert and Sommcrs 1995) Biermann 
has further suggested that the nearby radio galaxy M87 may be the source of the 
observed trans-GZK cosmic rays (see also Stecker 1968; Farrar and Piran 2000). Such 
an explanation would require one to invoke magnetic field configurations capable of 
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producing a quasi-isotropic distribution of > 10^° eV protons, making this hypothesis 
questionable. However, if the primary particles are nuclei, it is easier to explain a radio 
galaxy origin for the two highest energy events (Stecker and Salamon 1999; see section 
2.4). 

2.3.1. The Dead Quasar Origin Hypothesis It has been suggested that since all large 
galaxies are suspected to harbor supermassive black holes in their centers which may 
have once been quasars, fed by accretion disks which are now used up, that nearby 
quasar remnants may be the searched-for zcvatrons (Boldt and Ghosh 1999; Boldt and 
Lowenstein 2000). This scenario also has potential theoretical problems and needs to 
be explored further. In particular, it has been shown that black holes which are not 
accreting plasma cannot possess a large scale magnetic field with which to accelerate 
particles to relativistic energies (Ginzburg and Ozernoi 1964; Krolik 1999; Jones 2000). 
Observational evidence also indicates that the cores of weakly active galaxies have low 
magnetic fields (Falcke 2001 and references therein). Another proposed zevatron, the 
7-ray burst, is discussed in the next section. 

2.3.2. The Cosmological Gamma-Ray Burst Origin Hypothesis In 1995, it was 
hypothesized that cosmological 7-ray bursts (GRBs) could be the zevatron sources of 
the highest energy cosmic rays (Waxman 1995; Vietri 1995). It was suggested that 
if these objects emitted the same amount of energy in ultrahigh energy (~ 10"*^^ MeV) 
cosmic rays as in ~ MeV photons, there would be enough energy input of these particles 
into intergalactic space to account for the observed fiux. At that time, it was assumed 
that the GRBs were distributed uniformly, independent of redshift. 

Since 1997, X-ray, optical, and radio afterglows of more than two dozen GRBs 
have been detected leading to the subsequent identification of the host galaxies of these 
objects and consequently, their redshifts. To date, some 27 GRBs afterglows have been 
detected with a subsequent identification of their host galaxies. As of this writing, 26 
of the 27 are at moderate to high redshifts (> 0.36), with the highest one (GRB000131) 
lying at a redshift of 4.50. 

A good argument in favor of strong redshift evolution for the frequency of occurrence 
of the higher luminosity GRBs has been made by Mao and Mo (1998), based on the 
star-forming nature of the host galaxies. The host galaxies of GRBs appear to be sites 
of active star formation. The colors and morphological types of the host galaxies are 
indicative of ongoing star formation, as is the detection of Lya and [Oil] in several of 
these galaxies. Further evidence suggests that bursts themselves are directly associated 
with star forming regions within their host galaxies; their positions correspond to regions 
having significant hydrogen column densities with evidence of dust extinction. It now 
seems reasonable to assume that a more appropriate redshift distribution to take for 
GRBs is that of the average star formation rate. Results of the analyses of Schmidt 
(1999; 2001) also favor a GRB redshift distribution which follows the strong redshift 
evolution of the star formation rate. Thus, it now seems reasonable to assume that a 



more appropriate redshift distribution to take for GRBs is that of the average star 
formation rate, rather than a uniform distribution. If we thus assume a redshift 
distribution for the GRBs which follows the star formation rate, being significantly 
higher at higher redshifts, GRBs fail by at least an order of magnitude to account for 
the observed cosmic rays above 100 EeV (Stecker 2000). If one wishes to account for the 
GRBs above 10 EeV, this hypothesis fails by two to three orders of magnitude (Scully 
and Stecker 2002). Schmidt (2001) concludes that the local {z — 0) total energy release 
rate by all GRBs in the 7-ray range is < 10*^^ erg Mpc~^yr~^ whereas the required energy 
input rate in UHECRs above 10 EeV (10^^ cV) is 2.4 x 10^^ erg Mpc-^yr"^- § Even 
these numbers are most likely too optimistic, since they are based on the questionable 
assumption of the same amount of GRB energy being put into ultrahigh energy cosmic 
rays as in ~ MeV photons. 

2.3.3. Low Luminosity Gamma-Ray Bursts An unusual nearby Type Ic supernova, SN 
1998bw, at a redshift of 0.0085, has been identified as the source of a low luminosity 
burst, GRB980425, with an energy release which is orders of magnitude smaller than 
that for a typical cosmological GRB. Norris (2002) has given an analysis of the 
luminosities and space densities of such nearby low luminosity long-lag GRB sources 
which are identified with Type I supernovae. For these sources, he finds a rate per unit 
volume of 7.8 x 10^^ Mpc~^yr~^ and an average (isotropic) energy release per burst of 
1.3 xlO'^^ erg over the energy range from 10 to 1000 keV. The energy release per unit 
volume is then ~ 10^^ erg Mpc^'^yr"^. This rate is more than two orders of magnitude 
below the rate needed to account for the cosmic rays with energies above 10 EeV (see 
above) . 

2.4. The Heavy Nuclei Origin Scenario 

A more conservative hypothesis for explaining the trans-GZK events is that they were 
produced by heavy nuclei. Such nuclei would be required to come from astrophysical 
sources with a heavy composition and the conditions under which they were accelerated 
would have to preclude dissociation. Stecker and Salamon (1999) have shown that the 
energy loss time for nuclei starting out as Fe is longer than that for protons for energies 
up to a total energy of ~300 EeV (See Figure 6.) 

Stanev et al. (1995) and Biermann (1998) have examined the arrival directions of 
the highest energy events. They point out that the ~ 200 EeV event is within 10° of 
the direction of the strong radio galaxy NGC 315. This galaxy lies at a distance of 
only ~ 60 Mpc from us. For that distance, the results of Stecker and Salamon (1999) 
indicate that heavy nuclei would have a cutoff energy of ~ 130 EeV, which may be 

§ Recently, Vietri et al. (2003) have claimed as energy release rate of 1.1 x 10** erg Mpc~^yr~^. They 
obtain this rate by combining the total density given by Schmidt, which includes faint bursts, with 
a burst energy appropriate to the bright bursts. A more careful calculation from recent results gives 
2.5 X 10*^ erg Mpc'^^yr"^, a bit higher than the number quoted above, but still well below the required 
release rate for cosmic rays above 10 EeV. 
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Figure 6. Mean energy loss times for protons (Stecker 1968; Puget, Stecker and 
Bredelcamp 1976) and nuclei originating as Fe (Stecker and Salamon 1999). 

within the uncertainty in the energy determination for this event. The ~300 EeV event 
is within 12° of the direction of the strong radio galaxy 3C134. The distance to 3C134 
is unfortunately unknown because its location behind a dense molecular cloud in our 
own galaxy obscures the spectral lines required for a measurement of its redshift. 

An interesting new clue that we may indeed be seeing heavier nuclei above the 
proton-GZK cutoff comes from a recent analysis of inclined air showers above 10 EeV 
energy (Ave, et al. 2000). These new results favor proton primaries below the p-GZK 
cutoff energy but they appear to favor a heavier composition above the p-GZK cutoff 
energy. It will be interesting to see what future data from much more sensitive detectors 
will tell us. 

As can be seen from Figure 6, a spectrum of UHECR heavy nuclei should exhibit 
a "cutoff" at energies above 200 EeV (as opposed to 60 EeV for protons). Thus, an 
observed continuation of the UHECR spectrum to energies significantly above 200 eeV 
would rule out this origin hypothesis. 



2.5. Top-Down Scenarios: "Fraggers" 

A way to avoid the problems with finding plausible astrophysical zevatrons is to start at 
the top, i.e., the energy scale associated with grand unification, supersymmetric grand 
unification or its string theory equivalent. 

The modern scenario for the early history of the big bang takes account of the 
work of particle theorists to unify the forces of nature in the framework of Grand 
Unified Theories (GUTs) {e.g., , Georgi and Glashow 1974). This concept extends 
the very successful work of Nobel Laureates Glashow, Weinberg, and Salam in unifying 
the electromagnetic and weak nuclear forces of nature (Glashow 1960; Weinberg 1967; 
Salam 1968). As a consequence of this theory, the electromagnetic and weak forces would 
have been unified at a higher temperature phase in the early history of the universe and 
then would have been broken into separate forces through the mechanism of spontaneous 
symmetry breaking caused by vacuum fields which are known as Higgs fields. 

In GUTs, this same paradigm is used to infer that the electroweak force becomes 
unified with the strong nuclear force at very high energies of ~ 10^^ eV which occurred 
only ~ 10~^^ seconds after the big bang. The forces then became separated owing to 
interactions with the much heavier mass scale Higgs fields whose symmetry was broken 
spontaneously. The supersymmetric GUTs (or SUSY GUTs) provide an explanation 
for the vast difference between the two unification scales (known as the "Hierarchy 
Problem" ) and predict that the running coupling constants which describe the strength 
of the various forces become equal at the SUSY GUT scale of ~ 10^^ eV (Dimopoulos, 
Raby and Wilczek 1982). 

2.5.1. Topological Defects: Fossils of the Grand Unification Era Very heavy 
"topological defects" can be produced as a consequence of the GUT phase transition 
when the strong and electroweak forces became separated. These defects are locahzed 
regions of vacuum Higgs fields where extremely high densities of mass-energy are 
trapped. 

Topological defects in the vacuum of space are caused by misalignments of the 
heavy Higgs fields in regions which were causally disconnected in the early history of 
the universe. These are localized regions where extremely high densities of mass-energy 
are trapped. Such defects go by designations such as cosmic strings, monopoles, walls, 
necklaces (strings bounded by monopoles), and textures, depending on their geometrical 
and topological properties. Inside a topological defect vestiges of the early universe may 
be preserved to the present day. The general scenario for creating topological defects in 
the early universe was suggested by Kibble (1976). 

Superheavy particles or topological structures arising at the GUT energy scale 
M > 10^^ eV can decay or annihilate to produce "X-particles" (GUT scale Higgs 
particles, superheavy fermions, or leptoquark bosons of mass M.) In the case of strings 
this could involve mechanisms such as intersecting and intercommuting string segments 
and cusp evaporation. These X-particles will decay to produce QCD fragmentation 
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jets at ultrahigh energies, so I will refer to them as "fraggers" . QCD fraggers produce 
mainly pions, with a 3 to 10 per cent admixture of baryons, so that generally one can 
expect them to produce at least an order of magnitude more high energy 7-rays and 
neutrinos than protons. The same general scenario would hold for the decay of long- 
lived superheavy dark matter particles (see section 3.8), which would also be fraggers. It 
has also been suggested that the decay of ultraheavy particles from topological defects 
produced in SUSY-GUT models which can have an additional soft symmetry breaking 
scale at TeV energies ("flat SUSY theories") may help explain the observed 7-ray 
background flux at energies ~ 0.1 TeV (Bhattacharjee, Shafl and Stecker 1998). 

The number of variations and models for explaining the ultrahigh energy cosmic 
rays based on the GUT or SUSY GUT scheme (which have come to be called "top- 
down" models) has grown to be enormous and I will not attempt to list all of the 
numerous citations involved. Fortunately, Bhattacharjee and Sigl (2000) have published 
an extensive review with over 500 citations and I refer the reader to this review for 
further details of "top-down" models and references. The important thing to note here 
is that, if the implications of such models are borne out by future cosmic ray data, they 
may provide our flrst real evidence for GUTs. 

2.5.2. "Z-bursts" It has been suggested that ultra-ultrahigh energy (9(10 ZeV) 
neutrinos can produce ultrahigh energy fraggers by interactions with 1.9K thermal 
GBR neutrinos (Weiler 1982; Fargion et al. 1999; Weiler 1999), resulting in "Z-burst" 
fragmentation jets, again producing mostly pions. This will occur at the resonance 
energy E^es = 4[m,^(eV)]~^ ZeV. A typical Z boson will decay to produce ~2 nucleons, 
~20 7-rays and ~ 50 neutrinos, 2/3 of which are i/^'s. Gelmini and Kusenko (2000) 
have suggested a variant of model which involves both superheavy dark matter at high 
redshifts as an ultrahigh energy neutrino source (see Section 2.5.3) and subsequent Z- 
burst production at low redshifts. 

If the nucleons which are produced from Z-bursts originate within a few tens of 
Mpc of us they can reach us, even though the original ~ 10 ZeV neutrinos could have 
come from a much further distance. It has been suggested that this effect can be 
amplifled if our galaxy has a halo of neutrinos with a mass of tens of eV (Fargion, 
Mele and Sahs 1999; Weiler 1999). However, a neutrino mass large enough to be 
conflned to a galaxy size neutrino halo (Tremaine and Gunn 1979) or even a galaxy 
cluster size halo (Shafl and Stecker 1984) is now clearly ruled out by the results of the 
Wilkonson Microwave Anisitropy Probe {WMAP) (Spergel et al. 2003), combined with 
other cosmic microwave background data (Pearson et al. 2002; Kuo et al. 2002), data 
from the 2dF galaxy rcdshift survey (Colles et al. 2001; Elgaroy et al. 2002), together 
with the very small neutrino flavor mass differences indicated by the atmospheric and 
solar neutrino oscillation results (Fukuda et al. 1998, 1999; Smy 2002; Ahmad et al. 
2002; Bandyopadhyay et al. 2002; Gonzalez-Garcia and Nir 2003) imply that even the 
heaviest neutrino would have a mass in the sub-eV range, i.e., 0.03 eV < 777.3 < 0.24 eV 
(Bhattacharyya et al. 2003; Allen et al. 2003). The tritium decay spectral endpoint 
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limits on the mass of the electron neutrino (Weinheimer et al. 1999), are also consistent 
with this conclusion. This is much too small a mass for neutrinos to to be confined to 
halos of individual galaxies or evan galaxy clusters (Shafi and Stecker 1984). Because 
of the Pauli exclusion principle, the Fermi energy of such light neutrinos at the required 
number density to account for the halo dark matter would imply neutrino velocities far 
in excess of the escape velocity of these neutrinos from our galactic halo (Tremaine and 
Gunn 1979). 

The basic general problem with the Z-burst explanation for the trans-GZK events 
is that one needs to produce large fluxes of neutrinos with energies in excess of 10 
ZeV. If these are secondaries from pion production, this implies that the primary 
protons which produce them must have energies of hundreds of ZeV! Since we know 
of no astrophysical source which would have the potential of accelerating particles to 
energies even an order of magnitude lower (see section 4), a much more likely scenario 
for producing 10 ZeV neutrinos would be by a top-down process. The production rate of 
neutrinos from such processes is constrained by the fact that the related energy release 
into electromagnetic cascades which produce GeV range 7-rays is limited by the sateUite 
observations (see the review by Bhattacharjcc and Sigl 2000). This constraint, together 
with the low probability for Z-burst production, relegates the Z-burst phenomenon to a 
minor secondary role at best. 

2.5.3. Ultraheavy Dark Matter Particles: "Wimpzillas" The homogeneity and flatness 
of the present universe may imply that a period of very rapid expansion, called inflation, 
took place shortly after the big bang. The early inflationary phase in the history of the 
universe can lead to the production of ultrahigh energy neutrinos. The inflation of 
the early universe is postulated to be controlled by a putative vacuum fleld called the 
inflaton field. During inflation, the universe is cold but, when inflation is over, coherent 
oscillations of the inflaton fleld reheat it to a high temperature. While the inflaton field 
is oscillating, non-thermal production of very heavy particles ("wimpzillas") may take 
place. These heavy particles may survive to the present as a part of the dark matter. 
Their decays will produce ultrahigh energy particles and photons via fragmentation. 

It has been suggested that such particles may be the source of ultrahigh energy 
cosmic rays (Berezinsky et al. 1997; Kuz'min and Rubakov 1998; Blasi et al. 2002; 
Sarkar and Toldra 2002; Barbot and Drees 2002). The annihilation or decay of such 
particles in a dark matter halo of our galaxy would produce ultrahigh energy nucleons 
which would not be attenuated at trans-GZK energies owing to their proximity. A 
comparison of recent experimental constraints from dark matter nuclear recoil searches 
with predicted rates gives a lower limit on the mass of putative wimpzillas of 10^ EeV 
(Albuquerque amd Baudis 2003). 

A test of the halo wimpzilla hypothesis would be an arrival distibution which is 
skewed to favor the hemisphere in the direction of the inner galaxy and perhaps the inner 
galaxy itself. This would be an even larger effect in the case of wimpzilla annihilation 
(rather than decay), since the flux would then scale as the square of the wimpzilla 
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density (rather than hnearly). Since the galactic center is viewed from the southern 
hemisphere, the location of the AUGER detector will make it ideal for further testing 
the wimpzilla hypothesis. 

2.5.4- Halo Fraggers and the Missing Photon Problem Halo fragger models such as Z- 
burst and ultraheavy halo dark matter ("wimpzilla") decay or annihilation, as we have 
seen, will produce more ultrahigh energy photons than protons. These ultrahigh energy 
photons can reach the Earth from anywhere in a dark matter galactic halo, because, as 
shown in Figure 7, there is a "mini-window" for the transmission of ultrahigh energy 
cosmic rays between ~ 0.1 and ~ 10® EeV. 

Photon-induced giant air showers have an evolution profile which is significantly 
different from nucleon-induced showers because of the Landau-Pomeranchuk-Migdal 
(LPM) effect (Landau and Pomeranchuk 1953; Migdal 1956) and because of cascading in 
the Earth's magnetic field (Cillis et al. 1999) (see Figure 7). By taking this into account, 
Shinozaki, et al. (2002) have used the AGASA data to place upper limits on the photon 
composition of their UHECR showers. They find a photon content upper fimit of 28% 
for events above 10 EeV and 67% for events above 30 EeV at a 95% confidence level with 
no indication of photonic showers above 100 EeV. A recent reanalysis of the ultrahigh 
energy events observed at Haverah Park by Ave, et al. (2002) indicates that less than 
half of the events (at 95% confidence level) observed above 10 and 40 EeV are 7-ray 
initiated. An analysis of the highest energy Fly's Eye event {E = 300 EeV) shows it not 
to be of photonic origin as indicated in Figure 8. In addition, Shinozaki, et al. (2002) 
have found no indication of departures from isotropy as would be expected from halo 
fragger photonic showers, this admittedly with only 10 events in their sample. 

In order to solve the missing photon problem for halo fraggers, Chisholm and 
Kolb (2003) have suggested that a small violation of Lorentz invariance could allow 
ultrahigh energy photons to decay into electron-positron pairs, thus eliminating the 
photon component of the fraggcr-produced fiux. Photon decay by this mechanism 
was suggested by Coleman and Glashow (1999). The amount of Lorentz invariance 
required is within the observational limits obtained by Stecker and Glashow (2001) (see 
Section 3.10). However, the scenario suggested by Chisholm and Kolb, implies that 
neutrons would be the primary ultrahigh particles producing the giant air showers, 
again producing a halo anisotropy for which there is no present indication (Shinozaki, 
et al. 2002; KachelrieB and Semikoz 2003). 

2.6. Other New Physics Possibilities 

The GZK cutoff problem has stimulated theorists to look for possible solutions involving 
new physics. Some of these involve (A) a large increase in the neutrino-nucleon cross 
section at ultrahigh energies, (B) new particles, and (C) a small violation of Lorentz 
Invariance (LI). 
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Figure 7. The mean free path for ultrahigh energy 7-ray attenuation vs. energy. 
The curve for electron-positron pair production off the cosmic background radiation 
(CBR) is based on Gould and Schreder (1966). The two estimates for pair production 
off the extragalactic radio background are from Protheroe and Biermann (1996). The 
curve for double pair production is based on Brown, et al. (1973). The physics of pair 
production by single photons in magnetic fields is discussed by Erber (1966). This 
process eliminates all photons above ~ lO^** eV and produces a terrestrial anisotropy 
in the distribution of photon arrival directions above ~ lO^^eV. 
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Figure 8. The composite atmospheric shower profile of a 300 EeV photon- induced 
shower calculated with the Bethe-Heitler (dashed line) electromagnetic cross section 
and with the LPM effect taken into account (solid line, see text). The measured Fly's 
Eye profile, which fits the profile of a nucleonic primary, is shown by the data points 
(Halzen and Hooper 2002). 
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2.6.1. Increasing the Neutrino- Nucleon Cross Section at Ultrahigh Energies Since 
neutrinos can travel through the universe without interacting with the 2.7K CBR, 
it has been suggested that if the neutrino-nucleon cross section were to increase 
to hadronic values at ultrahigh energies, they could produce the giant air showers 
and account for the observations of showers above the proton-GZK cutoff. Several 
suggestions have been made for processes that can enhance the neutrino-nucleon cross 
section at ultrahigh energies. These suggestions include composite models of neutrinos 
(Domokos and Nussinov 1987; Domokos and Kovesi-Domokos 1988), scalar leptoquark 
resonance channels (Robinett 1988) and the exchange of dual gluons (Bordes, et al. 
1998).Burdman, Halzen and Gandhi (1998) have ruled out a fairly general class of these 
types of models, including those listed above, by pointing out that in order to increase 
the neutrino-nucleon cross section to hadronic values at ~ 10^° eV without violating 
unitarity bounds, the relevant scale of compositeness or particle exchange would have 
to be of the order of a GeV, and that such a scale is ruled out by accelerator experiments. 

More recently, the prospect of enhanced neutrino cross sections has been explored in 
the context of extra dimension models. Such models have been suggested by theorists 
to unify the forces of physics since the days of Kaluza (1921) and Klein (1926). In 
recent years, they have been invoked by string theorists and by other theorists as a 
possible way for accounting for the extraordinary weakness of the gravitational force, 
or, in other words, the extreme size of the Planck mass (Arkani-Hamed, Dimopoulos 
and Dvali 1999; Randall and Sundrum 1999). These models allow the virtual exchange 
of gravitons propagating in the bulk {i.e. in the space of full extra dimensions) while 
restricting the propagation of other particles to the familiar four dimensional space- 
time manifold. It has been suggested that in such models, (t(i/N) ~ [£',^/(10^'^eV)] 
mb (Nussinov and Schrock 1999; Jain, et al. 2000; see also Domokos and Kovesi- 
Domokos 1999). It should be noted that a cross section of ~ 100 mb would be necessary 
to approach obtaining consistency with the air shower profile data. Other scenarios 
involve the neutrino-initiated atmopheric production of black holes (Anchordoqui et 
al. 2002; Feng and Shapcrc 2002) and even higher dimensional extended objects, p- 
dimensional branes called "p-branes" (Ahn, Cavaglia and Olinto 2002; Anchordoqui, 
Feng and Goldberg 2002). Such interactions, in principle, can increase the neutrino 
total atmospheric interaction cross section by orders of magnitude above the standard 
model value. However, as discussed by Anchordoqui, Feng and Goldberg (2002), sub- 
mm gravity experiments and astrophysical constraints rule out total neutrino interaction 
cross sections as large as 100 mb as would be needed to fit the trans-GZK energy air 
shower profile data. Nonetheless, extra dimension models still may lead to significant 
increases in the neutrino cross section, resulting in moderately penetrating air showers. 
Such neutrino- induced showers should also be present at somewhat lower energies and 
provide an observational test for extra dimension TeV scale gravity models (Anchordoqui 
et al. 2001; Tyler, Ohnto and Sigl 2001). As of this writing, no such showers have been 
observed, putting an indirect constraint on fragger scenarios with TeV gravity models. 

There is also the possibility that the neutrino cross section extrapolated from lower 
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Figure 9. The air shower probability per incident Vt of energy 100 EeV as a function of 
neutrino cross section for upward moving air showers (UAS) and horizontal airshowers 
(HAS) assuming that the energy threshold for the detection of upward moving air 
showers (UAS) is 1 EeV (curve 1) and 10 EeV (curve 2). 

energies may be too high rather than too low (Dicus et al. 2001). Whatever the case, 
Kusenko and Weiler (2002) have argued that the ratio of upward moving air showers 
induced by ultrahigh energy z/^'s (from neutrino oscillations) having traveled through 
the Earth to horizontal neutrino induced air showers can be used to determine the 
neutrino- nucleon cross section at ultrahigh energies (see sections 4 and 5). The results 
of their calculation to determine this cross section are shown in Figure 9. 

2.6.2. New Particles The suggestion has also been made that new neutral hadrons 
containing a light gluino could be producing the trans-GZK events (Farrar 1996; Cheung, 
Farrar and Kolb 1998; Berezinsky, Kachelriefi and Ostapchenko 2002). While the 
invocation of such new particles is an intriguing idea, it seems unlikely that such particles 
of a few proton masses would be produced in copious enough quantities in astrophysical 
objects without being detected in terrestrial accelerators. There are strong accelerator 
constraints on light gluino production (Alavi-Harati, et al. 1999). 

One should note that while it is true that the GZK threshold for such particles 
would be higher than that for protons, such is also the case for the more prosaic heavy 
nuclei (see section 3.4). In addition, such neutral particles cannot be accelerated directly, 
but must be produced as secondary particles, making the energetics reqirements more 
difficult. 

2. 6. 3. Breaking Lorentz Invariance With the idea of spontaneous symmetry breaking 
in particle physics came the suggestion that Lorentz invariance (LI) might be weakly 
broken at high energies (Sato and Tati 1972). Although no real quantum theory of 
gravity exists, it was suggested that LI might be broken as a consequence of such a theory 
(Amelino-Camilia et al. 1998). A simpler formulation for breaking LI by a small first 
order perturbation in the electromagnetic Lagrangian which leads to a renormalizable 
treatment has been given by Coleman and Glashow (1999). Using this formalism, these 



authors have shown than only a very tiny amount of LI symmetry breaking is required to 
avoid the GZK effect by suppressing photomeson interactions between ultrahigh energy 
protons and the CBR. This LI breaking amounts to a difference of O{10~^^) between 
the maximum proton pion velocities. By comparison, Stecker and Glashow (2001) have 
placed an upper hmit of O{10~^^) on the difference between the velocities of the electron 
and photon, ten orders of magnitude higher than required to eliminate the GZK effect. 
Such a violation of Lorentz invariance might be produced by Planck scale effects (Aloisio, 
et al. 2000; Alfaro and Palma 2002, 2003). Lorentz invariance violation has also been 
invoked to solve the missing photon problem for halo fraggers (see Section 2.5.4). 

2.7. Is the GZK Effect All There Is? 

There is a remaining "dull" possibility. Perhaps the GZK effect is consistent with the 
data and is all there is at ultrahigh energies. The analysis of observations made with 
the HiRes monocular detector array shows only one event significantly above 100 EeV 
and appears to be consistent with the GZK effect (Abu-Zayyad et al. 2002; see Figure 
4 and section 3.7.) 

On the other hand. The strongest case for trans-GZK physics comes from the 
AGASA results. The AG ASA group, previously reported the detection of up to 17 
events with energy greater than or equal to ~ 100 EeV (Sasaki et al. 2001). They have 
recently reevaluated their energy determination and lowered this number. But they 
still claim a significant number of events at trans-GZK energies (Takeda, et al. 2003). 
The AGASA data indicate a deviation from pure GZK even if the number density of 
ultrahigh energy sources is weighted like the local galaxy distribution (Blanton, et al. 
2001). 

As mentioned earlier, De Marco et al. (2003) have argued that the discrepency 
between the AGASA and HiRes results is not statistically significant and that many 
more events are needed in order to determine whether or not there is a GZK effect. 
This will require data from future ground based detectors such as Auger and space 
based detectors such as EUSO and OWL (see Section 5). 

Further complicating the present observational situation, we note the fact that a 
fiuorescence detector such as HiRes, namely Fly's Eye, reported the highest energy 
event yet seen, it viz., E ~ 300 EeV. It is thus apparent that the experimental 
situation is interesting enough and the physics implications arc important enough to 
justify both more sensitive future detectors and the theoretical investigation of new 
physics and astrophysics. In this regard, we note that the Auger detector array will 
use both scintillators and fluorscence detectors (Zas 2001; see section 5.) Therefore, 
combined results from this detector array can help clarify the presently existing prima 
facie discrepency between the AGASA and HiRes results. 

It should also be noted that, even if the GZK effect is seen, top-down scenarios 
predict the reemergence of a new component at even higher energies (Aharonian, 
Bhattacharjee and Schramm 1992; Bhattacharjee and Sigl 2000). 
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2.8. Ultrahigh Energy Event Signatures 

Future data which will be obtained with new detector arrays and satellites (see next 
section) will give us more clues relating to the origin of the trans-GZK events by 
distinguishing between the various hypotheses which have been proposed. 

A zevatron origin ("bottom-up" scenario) will produce air-showers primarily from 
primaries which are protons or heavier nuclei, with a much smaller number of neutrino- 
induced showers. The neutrinos will be secondaries from the photomeson interactions 
which produce the GZK effect (Stecker 1973; 1979; Engel, Seckel and Stanev 2001 and 
references therein). In addition, zevatron events may cluster near the direction of the 
sources. 

A "top-down" (GUT) origin mechanism will not produce any heavier nuclei and 
will produce more ultrahigh energy neutrinos than protons. This was suggested as a 
signature of top-down models by Aharonian, Bhatacharjee and Schramm (1992). Thus, 
it will be important to look for the neutrino-induced air showers which are expected 
to originate much more deeply in the atmosphere than proton-induced air showers and 
are therefore expected to be mostly horizontal showers. Looking for these events can 
most easily be done with a satellite array which scans the atmosphere from above (see 
Section 4.) 

Top-down models also produce more photons than protons. However, the mean 
free path of these photons against pair-production interactions with extragalactic low 
frequency radio photons from radio galaxies is only a few Mpc at most (Protheroe and 
Biermann 1996). The subsequent electromagnetic cascade and synchrotoron emission of 
the high energy electrons produced in the cascade dumps the energy of these particles 
into much lower energy photons (Wdowczyk, Tkaczyk and Wolfcndale 1972; Stecker 
1973). However, the photon-proton ratio is an effective tool for testing halo fragger 
models (See sect. 3.9.) 

Another characteristic which can be used to distinguish between the bottom-up 
and top-down models is that the latter will produce much harder spectra. If differential 
cosmic ray spectra are parametrized to be of the form F oc E~^, then for top-down 
models F < 2, whereas for bottom-up models F > 2. Also, because of the hard source 
spectrum in the "top-down" models, they should exhibit both a GZK suppression and 
a pileup just before the GZK energy. 

If Lorentz invariance breaking is the explanation for the missing GZK effect, the 
actual absence of photomeson interactions should result the absence of a pileup effect 
as well. 

3. The Highest Energy Gamma Rays 

The highest energy 7-rays observed to date were produced by the Vela pulsar and by 
PSR 1706-44, by supernova remnants in our galaxy and by extragalactic sources known 
as blazars. Blazars are a class active galaxies believed to have supermassive black 
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holes in their nuclei which gravitationally power jets which produce massive amounts 
of nonthermal radiation. They are distinguished by the condition that their jets are 
pointed almost directly at us, producing interesting relativistic effects such as rapid 
time variability and Lorentz boosted radiation and also superluminal motion in many 
cases (Jorstad et al. 2001). 

3.1. Pulsars 

There are two theoretical models which have been proposed to account for the origin of 
pulsed 7-ray emission in pulsars, viz., "the outer gap model" , where particle acceleration 
occurs in the outer magnetosphere of the pulsar (Cheng, Ho and Ruderman 1986), and 
"the polar cap model" , where particle acceleration occurs near the magnetic polar cap 
of the pulsar (Daugherty and Harding 1996). 

In the outer gap model, in the case of the Crab pulsar, the resulting electron- 
positron pairs can generate TeV 7-rays through the synchrotron-self-Compton (SSC) 
mechanism. In SSC emission models, the source has a natural two-peaked spectral 
energy distribution. The lower energy peak is produced by synchrotron radiation of 
relativistic electrons accelerated in the source and the higher energy peak is produced by 
these same electrons upscattering the synchrotron component photons to TeV energies 
by Compton interactions. In the case of other pulsars, the low energy photons may 
come from curvature radiation of electrons and positrons as they follow the magentic 
field lines of the pulsar. This curvature radiation would then be Compton upscattered. 

In the polar cap model, the strong magnetic field near the pulsar would cut off 
any TeV photons emitted near the surface of the pulsar via single-photon electron- 
positron pair production off the magnetic field (Erbcr 1966). Thus, the detection of 
TeV pulsed emission from ordinary pulsars would favor the outer gap model. However, 
for miUisecond pulsars, whose 5-fields are 3-4 orders of magnitude lower than is the 
case for regular pulsars, the high energy cutoff from single-photon pair production is 3-4 
orders of magnitude higher in energy (Bulik et al. 2000) .Thus, TeV emission may be 
detected in these sources in the future even in the case of the polar cap 7-ray production. 

3.2. Supernova Remnants 

Very high energy (TeV) 7-rays have been reported from several supernova remnants, 
viz., the Crab Nebula (Weekes et al. 1989), the Vela pulsar wind nebula (Yoshikoshi et al. 
1997), and the nebulae associated with PSR 1706-44 (Kifune et al. 1995), Cassiopeia A 
(Piihlhofer et al. 1999), SN1006 (Tanimori et al. 1989), and RXJ1713.7-3946 (Muraishi 
et al. 2000). 

The SSC mechanism has been invoked to account for the TeV emission in the Crab 
Nebula (de Jager and Harding 1992). Figure 10 from this paper shows the observed 7-ray 
spectrum of the Crab Nebula together with the theoretical fits for the synchrotron and 
SSC Compton components. In this interpretation, the synchrotron component extends 
to an energy of 0.1 GeV, with the Compton component dominating at higher energies. 
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Figure 10. The observed 7-ray spectrum of the Crab Nebula with synchrotron and 
Compton components shown. 



Another mechanism which has been proposed to explain the TeV emission in 
supernova remnants is the Compton scattering of very high energy relativistic electrons 
off the 2.7 K cosmic background radiation. This has been proposed for specifically for 
the remnant SN1006 {e.g., , Pohl et al. 1996; Mastichiadas and de Jagcr 1996). 

There is also the important possibility that TeV 7-rays could be produced 
in supernova remnants by accelerated highly relativistic protons interacting with 
interstellar gas nuclei in the vicinity of the remnant and producing very high energy 7- 
rays through the mechanism of the production and decay of neutral pions (tt^'s) (Drury, 
Aharonian and Volk 1994; Gaisser, Protheroe and Stanev 1998). Almost 70 years ago, 
Baade and Zwicky (1934) first proposed that supernovae could provide the energy for 
accelerating cosmic rays in our galaxy. This hypothesis gained observational support two 
decades later when Shklovskii (1953) proposed that relativistic electrons radiating in the 
magnetic field of the Crab Nebula produced its optical continuum radiation. Indirect 
support for proton acceleration in supernovae came from 7-ray observations in the 1970s 
(Stecker 1975; 1976). It now appears that evidence for the hadronic production 
mechanism may have been found for the source RXJ1713. 7-3946 (Enomoto et al. 2002). 
While the evidence is not definitive (Butt et al. 2002), we can hope for a resolution 
with future observations. The detection of high energy neutrinos from RXJ1713. 7-3946 
would serve as a smoking gun for hadronic producion since these neutrinos would be 
produced by the decay of tt^ mesons (Alvarez-Mianiz and Halzen 2002). 

3.3. Blazars 

Blazars were first discovered as the dominant class of extragalactic high energy 7-ray 
sources by the EGRET detector on the Compton Gamma Ray Observatory (CGRO). 
Because they are at large extragalactic distances, their spectra are predicted to be 
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modified by strongly redshift dependent absorption effects caused by interactions of 
tfiese 7-rays witli photons of tlie intergalactic IR-UV background radiation (Stecker, de 
Jager and Salamon 1992). 

Tlie liigliest energy extragalactic 7-ray sources are tliose blazars known as X-ray 
selected BL Lac objects (XBLs), or alternatively as high frequency BL Lac objects 
(HBLs). They are expected to emit photons in the multi-TeV energy range (Stecker, de 
Jager and Salamon 1996), but only the nearest ones are expected to be observable, the 
others being hidden by intergalactic absorption (Stecker et al. 1992). 

There are now ~70 "grazars" (7-ray blazars) which have been detected by the 
EGRET team at GeV energies (Hartman et al. 1999) These sources, optically violent 
variable quasars and BL Lac objects, have been detected out to a redshift of ~2.3. 

In addition, several TeV grazars have been discovered at low redshifts [z < 0.13), 
viz., Mkn 421 at z = 0.031 (Punch et al. 1992), Mkn 501 at z = 0.034 (Quinn et al. 
(1996), 1ES2344+514 at z = 0.44 (Catanese et al. 1998), PKS 2155-304 at z = 0.117 
(Chadwick et al. 1998), and 1ES1426+428 (a.k.a. H1426+428) at z = 0.129 (Aharonian 
et al. 2003; Horan et al. 2002). These sources all fit the two candidate source criteria 
of Stecker et al. (1996) of closeness and a high frequency synchrotron peak which is 
prominant in X-ray emission. Recently, another criterion has been suggested for TeV 
candidate sources, viz. a relatively large flux in the synchrotron peak (Costamante 
and Ghisellini 2001). The closeness criterion has to do with extragalactic TeV 7- 
ray absorption by pair production (Stecker et al. 1992). The other two criteria have 
to do with the synchrotron-self- Compton (SSC) mechanism believed to be primarily 
responsible for the TeV emission of the HBL sources. In the SSC models, an HBL 
blazar has a two-peaked spectral energy distribution (see pulsar section above) with the 
lower energy peak in the radio to X-ray range and the higher energy peak in the X-ray 
to multi-TeV 7-ray range. 

Very high energy 7-ray beams from blazars can be used to measure the intergalactic 
infrared radiation field, since pair-production interactions of 7-rays with intergalactic 
IR photons will attenuate the high-energy ends of blazar spectra (Stecker et al. 1992). 
In recent years, this concept has been used successfully to place upper limits on the the 
diffuse infrared radiation background (DIRB) (Stecker and de Jager 1993, 1997; Dwek 
1994; Stanev and Pranceschini 1997; Biller et al. 1998). Determining the DIRB, in turn, 
allows us to model the evolution of the galaxies which produce it. As energy thresholds 
are lowered in both existing and planned ground-based air Cerenkov light detectors and 
with the launch of the GLAST (Gamma-Ray Large Area Space Telescope) in 2006, 
cutoffs in the 7-ray spectra of more distant blazars are expected, owing to extinction by 
the DIRB. These can be used to explore the redshift dependence of the DIRB (Salamon 
and Stecker 1998 (SS98)). 

In addition, as we shall discuss, the very existence of an absorption feature in the 
spectra of TeV blazars provides a sensitive test for the exactness of Lorentz invariance. 
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3.4- The Diffuse Low Energy Photon Background and Extragalactic Gamma Ray 
Absorption 

The formulae relevant to absorption calculations involving pair-production interactions 
with redshift factors are given in Stecker et al. (1992). For 7-rays in the TeV energy 
range, the pair-production cross section is maximized when the soft photon energy is in 
the infrared range: 

\{E^) ~ Ae^^ = l.'2^E^,TeV l^m (1) 

where \e = h/ (mec) is the Compton wavelength of the electron. For a 1 TeV 7-ray, this 
corresponds to a soft photon having a wavelength ~ 1 /xm. (Pair-production interactions 
actually take place with photons over a range of wavelengths around the optimal value 
as determined by the energy dependence of the cross section.) If the emission spectrum 
of an extragalactic source extends beyond 20 TeV, then the extragalactic infrared field 
should cut off the observed spectrum between ~ 20 GeV and ~ 20 TeV, depending on 
the redshift of the source (Stecker and Salamon 1997; SS98). 

Several attempts have been made to infer the IR SED (spectral energy distribution) 
of the DIRB either from model calculations or observations. (See Hauser and Dwek 
(2001) for the latest review.) Such information can be used to calculate the optical 
depth for TeV range photons as a function of energy and redshift. Figure 11 summarizes 
the observationally derived values for the extragalactic optical-UV, near-IR and far- 
IR fluxes which now exist. We will refer to the totahty of these fluxes as the 
extragalactic background light (EBL) Unfortunately, foreground emission prevents the 
direct detection of the EBL in the mid-IR wavelength range. (Sec discussion in Hauser 
and Dwek 2001.) However, other theoretical models such as those of Tan, Silk and 
Balland (1999), Rowan- Robinson (2000) and Xu (2000) predict fairly flat SEDs in the 
mid-infrared range with average flux levels in the 3 to 4 nW m"^sr as do the Malkan 
and Stecker (2001) models shown in Figure 11. These flux levels are also consistent with 
the indirect mid-IR constraints indicated by the box in Figure 11. (These constraints 
are summarized by Stecker (2001).) || 

The two Malkan and Stecker (2001) SED curves for the DIRB, extended into the 
optical-UV range by the hybrid model of de Jager and Stecker (2002) (DS02), give a 
reasonable representation of the EBL in the UV to far-IR. Other EBL models in the 
literature whose flux levels roughly fit the present data and have the same spectral 
characteristics, i.e., a stellar optical peak, a far-IR dust emission peak, and a mid-IR 
valley which allows for some warm dust emission (see review of Hauser and Dwek 2001), 
should give similar results on the optical depth of the near Universe t{E,z) to high 
energy 7-rays. 

II We note that the COBE-DIRBE group has argued that a real flux derived from the COBE-DIRBE 
data at 100 /iin, as claimed by Lagache et al. (2000), is untenable because isotropy in the residuals 
(after foreground subtractions) could not be proven. Dwek et al. (1998) have concluded that only a 
conservative lower limit of 5 m~^sr~^ could be inferred at 100 ^,m. 
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Figure 11. The SED of the EBL (see text for references and descriptions). All 
error bars given at the ±2a level. The Malkan and Stecker (2001) fast evolution 
model is shown by the upper, thick solid curve between log^g ^/ = 11.8 to 13.8) and 
their baseline model is shown by the thick dashed line over the same frequency range. 
Convergent Hubble Deep Field galaxy counts (Madau and Pozzetti 2001): open circles; 
Ground-based galaxy counts limits at 2.2/im (see text): thick vertical bar marked 
"2.2/i"; COBE-DIRBE photometric sky residuals (Wright and Reese 2000): small 
open squares; TeV 7-ray-based upper limits: thick box marked "TeV" (see text for 
references); ISOCAM lower limit at 15/im: upward arrow marked 15/x (Elbaz et al. 
1999); COBE-DIRBE far-IR sky residuals (Hauser et al. 1998): large open squares 
at 140 /im and 240 yum; COBE-DIRBE data recalibrated using the COBE- FIRAS 
calibration (see text): large solid squares without error bars; ISOCAM 170 jim. flux 
(Kiss et al. 2001): solid circle; COBE-FIRAS far-IR sky residuals (Fixsen, et al. 
1998): thick dot-dash curve with ^ 95% confidence band (thin dot-dash band); lOO/xm 
COBE-DIRBE point (Lagache et al. 2000): small solid square; flux at 3.5 /xm from 
COBE-DIRBE (Dwek and Arendt 1998): sohd diamond. 

DS02 rcdcrivcd the optical depth of the universe to high energy 7-rays as a function 
of energy and redshift for energies betweeen 50 GeV and 100 TeV and redshifts between 
0.03 and 0.3 using their derived hybrid EBL SEDs. Figure 12 shows t{E, z) calculated 
using the basehne and fast evolution hybrid models of DS02 for 7-ray energies down to 
~ 50 GeV, which is the approximate threshold energy for meaningful image analyses in 
next generation ground based 7-ray telescopes such as MAGIC, H.E.S.S. and VERITAS. 
Where they overlap, the DS02 results are in good agreement with the metallicity 
corrected results of SS98, which give the optical depth of the universe to 7-rays out 
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Figure 12. The optical depth for 7-rays above 50 GeV given for redshifts between 
0.03 and 0.3 (as labelled) calculated using the medium (dashed lines) and fast evolution 
(solid hues) SED of Malkan and Stecker (2001). 

to a redshift of 3 and extend to lower energies (See Section 2.7.)^ 
3.5. The Nearby TeV Blazar Source Mkn 501 

As can be seen in Figure 12, intergalactic space is predicted to become opaque for 7-rays 
having an energy greater than ~ 3 TeV originating at a redshift of ~ 0.03, as is the case 
with the TeV 7-ray sources MknSOl and Mkn 421. 

The nearby blazar Mkn 501 is of particular interest because a very detailed 
determination of its spectrum was obtained by observing it while it was strongly flaring 
in 1997. The spectrum observed at that time by the HEGRA air Cerenkov telescope 
system (Aharonian, et al. 2001a) extended to energies greater than 20 TeV, the highest 
energies yet observed from an extragalactic source. 

Using this observational data, DS02 derived the intrinsic 7-ray spectrum of Mkn 501 
during its 1997 flaring state by compensating for the effect of intergalactic absorption. 

% Owing to a calculational inaccuracy, the curves shown for z = 0.03 are too high for the assumed 
SEDs by roughly 30%. This has been corrected in Konopclko et al. (2003) thanks to it being pointed 
out by Eli Dwek. However, given that the uncertainties in the SEDs are of this order and given that 
the actual redshift of Mkn 501 is 0.034, this does not represent a significant difference. 
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Figure 13. The observed spectrum and derived intrinsic spectrum of Mkn 501. 

The observed spectral data are as measured by HEGRA (sohd circles) and Whipple 
(solid squares). The upper points are the absorption corrected data points (marked 
"INTRINSIC") using our fast evolution hybrid EBL (upper data set and solid curve 
fit) and baseline hybrid EBL (lower data set with dashed curve fit). 

They found that the time averaged spectral energy distribution of Mkn 501 while flaring 
had a broad, flat peak in the ~5-10 TeV range which corresponds to the broad, flat time 
averaged X-ray peak in the ~ 50-100 keV range observed during the flare. The spectral 
index of the derived intrinsic differential photon spectrum for Mkn 501 at energies below 
~2 TeV was found to be ~1.6-1.7. This corresponds to a time averaged spectral index 
of 1.76 found in soft X-rays at energies below the X-ray (synchrotron) peak (Petry et al. 
2000). These results appear to favor an SSC origin for the TeV emission together with 
jet parameters which are consistent with time variability constraints within the context 
of a simple SSC model. The similarity of the soft X-ray and ~TeV spectral indices 
of the two components of the source spectrum implies that 7-rays below ~ 1 TeV are 
produced in the Thomson regime by scattering off synchrotron photons with energies in 
the optical-IR range. On the other hand, 7-rays near the ~ 7 ± 2 TeV Compton peak 
appear to be the result of scattering in the Klein-Nishina range. 

The observed spectrum between 0.56 and 21 TeV was obtained from 
contemporaneous observations of Mkn 501 by the HEGRA group (Aharonian et al. 
2001a) and the Whipple giouY> (Krennrich et al. 1999). These observations are consistent 
with each other (within errors) in the overlapping energy range between 0.5 and 10 
TeV, resulting in a single spectrum extending over two decades of energy, marked 
"OBSERVED" in Figure 13. 
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Using both the Whipple and HEGRA data and correcting for absorption by 
multiplying by e'^^^^ evaluated at 2; = 0.034 with their newly derived values for the 
opacity, DS02 derived the intrinsic spectrum of Mkn 501 over two decades of energy. 
This is given by the data points and two curves marked "INTRINSIC" in Figure 13. The 
upper of these curves corresponds to the fast evolution case; the lower curve corresponds 
to the baseline model case. 

Fossati et al. (2000) has suggested a parameterization to describe smoothly curving 
blazar spectra. This parameterization is of the form 



A spectrum of this form changes gradually from a spectral index of Fi to an index of F2 
when the energy E increases through the break energy Eb- The parameter / describes 
the rapidity ( "fastness" ) of the change in spectral index over energy. DS02 apphed the 
formalism of Fossati et al. (2000) to their Mkn 501 intrinsic flare spectrum and found 
best-fits for the parameters K, Eb, f, Fi and F2 after correcting the observed spectrum 
for intergalactic absorption. 

Whereas the low energy spectral index Fi was found to be well constrained, 
the higher energy index F2 is unconstrained. The SED peaks at Em ~ 8 — 9 TeV 
(independent of the unconstrained F2) in the case where the fast evolution EBL is 
assumed and that Em ~ 5 TeV if the basehne EBL is assumed. (See Figure 13) 

Since Mkn 501 is a giant elliptical galaxy with httle dust, it is reasonable to assume 
that the galaxy itself does not produce enough infrared radiation to provide a significant 
opacity to high energy 7-rays. Such BL Lac objects have little gas (and therefore most 
likely little dust) in their nuclear regions. It also appears that 7-ray emission in blazars 
takes place at super luminal knots in the jet downstream of the core and at any putative 
accretion disk. So, if the EBL SED of DS02 is approximately correct, it is reasonable to 
assume that the dominant absorption process is intergalactic and that pair-production 
in the jet in negligible. This hypothesis is supported by the fact that the high energy 
7-ray SED did not steepen during the flare. This implies that the optical depth given 
by eq. (4) is less than unity out to the highest observed energy £^ ~ 20 TeV. Thus, it 
appears that the high energy turnover in the observed Mkn 501 7-ray spectrum above 
~ 10 TeV can be understood solely as a result of intergalactic absorption. 

The intrinsic SED of Mkn 501 derived by DS02 is quite flat in the multi-TcV range as 
shown in Figure 3. This is in marked contrast to the dramatic turnover in its observed 
SED. This is strong evidence that the observed spectrum shows just the absorption 
effect predicted. We will see that the spectral observations of other blazars, although 
not nearly as good in most cases, also exhibit evidence of intergalactic absorption. 

3.6. Absorption in the Spectra of other Blazars 

Observations of Mkn 421, the closest TeV blazar which has a redshift very similar to 
that of Mkn 501, were made by the Whipple group up to an energy of 17 TeV (Krennrich 
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Figure 14. Predicted differential absorbed spectrum, for PKS 2155-304 (solid line) 
assuming an differential source spectrum (dashed line) normalized to the integral 
flux given by Chadwick et al. (1998). 

et al. 2002). Their results indicate a turnover in the spectrum at energies above ~4 TeV 
caused by extragalactic absorption, quite similar to that observed in the spectrum of 
Mkn 501 (see above). 

As to sources at somewhat higher rcdshifts, Stecker (1999) considered the blazar 
source PKS 2155-304, located at a moderate redshift of 0.117, which has been reported 
by the Durham group to have a flux above 0.3 TeV of ~ 4 x 10~^^ cm~^ s~^, close 
to that predicted by a simple SSC model. Using absorption results obtained by Stecker 
and de Jager (1998) and assuming an E~'^ source spectrum, Stecker (1999) predicted 
an absorbed (observed) spectrum as shown in Figure 14. As indicated in the figure, it 
was predicted that this source should have its spectrum steepened by ~ 1 in its spectral 
index between ~ 0.3 and ~ 3 TeV and should show a pronounced absorption turnover 
above ~ 6 TcV. 

There arc recent observations of the spectrum of another TcV blazar lES 1426+428 
at a redshift of 0.129, not too different from that of PKS 2155-304. Interestingly, it 
has been found by three groups that the spectrum of this source is quite steep, with 
a photon differential spectral index greater than 3 (Aharonian et al. 2002; Petry et al. 
2002; Djannati-Atai, et al. 2002). This steep spectrum, similar to that predicted by 
Stecker (1999) for PKS2155-304 at a redshift of 0.117, is more evidence for extragalactic 
absorption. However, the spectral data for this source are not as high a quality as for 
the more nearby sources Mkn 421 and Mkn 501 so that the details of the absorption 
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cannot be accurately determined. 

3.1. The Gamma-Ray Opacity at High Redshifts 

Salamon and Stecker (1998) (SS98) have calculated the 7-ray opacity as a function of 
both energy and rcdshift for redshifts as high as 3 by taking account of the evolution 
of both the SED and emissivity of galaxies with redshift. In order to accomplish this, 
they adopted the recent analysis of Fall et al. (1996) and also included the effects 
of metallicity evolution on galactic SEDs. They then gave predicted 7-ray spectra 
for selected blazars and extend our calculations of the extragalactic 7-ray background 
from blazars to an energy of 500 GeV with absorption effects included. Their results 
indicate that the extragalactic 7-ray background spectrum from blazars should steepen 
significantly above 15-20 GeV, owing to extragalactic absorption. In addition, there 
are most likely cutoffs in the intrinsic spectra of the most powerful 7-ray blazars above 
this energy range. Thus, future observations of such a steepening in the background 
spectrum by GLAST would provide evidence for the blazar origin hypothesis for the 
7-ray background radiation. 

SS98 calculated stellar emissivity as a function of redshift at 0.28 /xm, 0.44 /xm, and 
1.00 /xm, both with and without a metallicity correction. Their results agree well with 
the emissivity obtained by the Canada-Prance Redshift Survey (Lilly et al. 1996) over 
the redshift range of the observations {z < 1). The stellar emissivity in the universe 
is found to peak at 1 < 2; < 2, dropping off steeply at lower reshifts and is roughly 
constant higher redshifts {e.g., Steidel 1999). Indeed, Madau and SchuU (1996) have 
used observational data from the Hubble Deep Field to show that metal production 
has a similar redshift distribution, such production being a direct measure of the star 
formation rate. (See also Pettini et al. 1994.) 

The optical depth of the universe to 7-rays as a function of energy at various 
redshifts out to 2; = 3 which was derived by SS98 is shown in Figure 15. 

3.8. The Effect of Absorption on the Spectra of Blazars and the Gamma-Ray 
Background 

With the 7-ray opacity t{Eq,z) calculated out to z — ?> (see previous section), the 
cutoffs in blazar 7-ray spectra caused by extragalactic pair production interactions with 
stellar photons can be predicted. The left graph in Figure 16 from SS98 shows the effect 
of the intergalactic radiation background on a few of the grazars observed by EGRET, 
viz., 1633+382, 3C279, 3C273, and Mkn 421, assuming that the mean spectral indices 
obtained for these sources by EGRET extrapolate out to higher energies attenuated 
only by intergalactic absorption. In considering this figure, it should be noted that 
observed cutoffs in many grazar spectra may also be affected by natural cutoffs in their 
source spectra (Stecker, de Jager and Salamon 1996) and intrinsic absorption may also 
be important in some sources (Protheroe and Biermann 1996). 
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Figure 15. The optical depth for 7-rays as a function of energy for various redshifts 
calculated using the metallicity correction of SS98. 



Figure 17 shows the background spectrum predicted from unresolved blazars 
(Stecker and Salamon 1996a; SS98) , compared with the EGRET data (Sreekumar et 
al. 1998). Note that the predicted spectrum steepens above 20 GeV, owing to 
extragalactic absorption by pair-production interactions with radiation from external 
galaxies, particularly at high redshifts. 

Again we note that the predicted background above 10 GeV from unresolved blazars 
is uncertain because many blazars are expected to have intrinsic cutoffs in their 7- 
ray production spectra and by intrinsic 7-ray absorption within such sources is also a 
possibility. Thus, above 10 GeV the calculated background flux from unresolved blazars 
shown in Figure 17 may actually be an upper limit. Whether cutoffs in grazar spectra 
are primarily caused by intergalactic absorption can be determined by observing whether 
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Figure 16. The effect of intergalactic absorption by pair-production on the power- 
law spectra of four prominent grazars: 1633+382 (z = 1.81), 3C279 {z = 0.54), 3C273 
(z ~ 0.15), and Mkn 421 (z = 0.031) assuming an extrapolation of the spectral indeces 
for these sources measured by EGRET to high energy. 

the grazar cutoff energies have the type of redshift dependence predicted here. 
3.9. Constaints on the Redshifts of Gamma-Ray Bursts 

The results of the SS98 absorption calculations can also be used to place limits on the 
redshifts (or distances) of 7-ray bursts. On 17 February 1994, the EGRET telescope 
observed a 7-ray burst which contained a photon of energy ~ 20 GeV (Hurley et al. 
1994). If one adopts the opacity results which include the metallicity correction, the 
highest energy photon in this burst would be constrained probably to have originated 
at a redshift less than ~2. 

On 17 April 1997, the ground based "Milagrito" detector observed the burst 
GRB970417a with an effective threshold of energy of 0.65 TeV (Atkins et al. 2003) .Using 
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Figure 17. The extragalactic 7-ray background spectrum predicted by the unresolved 
blazar model of Stecker and Salamon (1996) with absorption included, calculated for a 
mean EGRET point-source sensitivity of 10~^ cm~^s~^, compared with the EGRET 
data on the 7-ray background (Sreekumar et al. 1998). The solid (dashed) curves are 
calculated with (without) the metallicity correction function (from SS98). 

the opacity curves derived by DS02 as shown in Figure 12, we can obtain constraints 
on the maximum redshift of this GRB. Since the observed burst contained photons of 
energies down to 0.65 TeV, the maximum redshift of GRB970417a was z^ax — 0.1. 

Future detectors such as GLAST (Bloom 1996; Michelson 2001) should be able to 
place redshift constraints on bursts observed at higher energies. Such constraints may 
further help to identify the host galaxies of 7-ray bursts. 

3.10. Constraints on the Violation of Lorentz Invariance. 

Lorentz invariance violation can be described quite simply in terms of different maximal 
attainable velocities of different particle species as measured in the preferred frame 



(Coleman and Glashow 1999). Following the well-defined formalism for LI violation 
discussed by Coleman and Glashow (1999), (see also Colladay and Kostelecky 1998), 
the maximum attainable velocity of an electron need not equal the in vacua velocity of 
light, i.e., Ce ^ Cj. The physical consequences of this violation of LI depend on the sign 
of the difference. Defining 



Stecker and Glashow (2001) consider the two cases of positive and negative values of S 
separately. 

Case I: If Ce < {S < 0), the decay of a photon into an electron-positron pair is 
kinematically allowed for photons with energies exceeding 



The decay would take place rapidly, so that photons with energies exceeding ii^max could 
not be observed either in the laboratory or as cosmic rays. Prom the fact that photons 
have been observed with energies > 50 TeV from the Crab nebula (Tanimori et al. 
1998), it follows from eq.(9) that E^^ > 50 TeV, or that -5 <2x 10'^^. 

Case II: Here we are concerned with the remaining possibility, where Cg > Cj 
{S > 0) and electrons become superluminal if their energies exceed £^max/ V^- Electrons 
traveling faster than light will emit light at all frequencies by a process of 'vacuum 
Cerenkov radiation.' This process occms rapidly, so that superluminal electron energies 
quickly approach -Emax/ V^- However, because electrons have been seen in the cosmic 
radiation with energies up to ~ 1 TeV (Nishimora et al. 1980), one obtains an upper 
hmit on S in this case of 3 x 10~^^. This hmit is two orders of magnitude weaker than 
the limit obtained for Case I. However, if the observed ~TeV 7-ray emission from the 
Crab Nebula is produced by very high energy electrons via the SSC mechanism (de 
Jager and Harding 1992), then electrons of energy > 50 TeV are required to produce 
the observed 50 TeV 7-rays in the nebula. Assuming that this is the case, one obtains 
a more indirect upper limit on 5 of 1 x 10~^^. 

Stecker and Glashow (2001) have also shown how stronger bounds on 6 can be set 
using observations of very high energy cosmic ray photons. For case I, the discussion 
is trivial: The mere detection of cosmic 7-rays with energies greater that 50 TeV from 
sources in our galaxy would improve the bound on S. For case II, if LI is broken so 
that Ce > c-y, the threshold energy for the pair production process 7 -|- 7 — > e"*" -|- e~ is 
altered.The square of the four-momentum becomes 



where e is the energy of the low energy (infrared) photon and 6 is the angle between 
the two photons. The second term on the left-hand-side comes from the fact that 



Ce = c^{l + S) , < |5| < 1 , 
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2eE^(l - cos^) - 2E^S = A-f^ml > Ami 



(5) 



For head-on collisions (cos^ — —1) the minimum low energy photon energy for pair 
production becomes 

emin^ml/E^ + {E^5)/2 (6) 

It follows that the condition for a significant increase in the energy threshold for pair 
production is 

Er-f > ^'max or equivalently S > 2ml/ ■ i"^) 

Extragalactic photons with the highest energies yet observed originated in a 
powerful flare coming from Mkn 501 (Aharonian et al. (1999) .We have seen that the 
Mkn 501 observations indicate that its multi-TeV spectrum is consistent with what one 
would expect from intergalactic absorption (SD02). Since there is no significant decrease 
in the optical depth of the universe at the distance of Mkn 501 for E^ < 20 TeV, it then 
follows from eq. (5) that 6 < 2{me/ E^y = 1.3 x 10~^^. This constraint is two orders of 
magnitude stronger than that obtained from the direct cosmic-ray electron data. 

Future detection of galactic 7-rays with energies greater than 50 TeV would 
strengthen the bound on 6 for Case I. For Case II, the detection of cosmic 7-rays 
above 100(1 -I- Zg)'"^ TeV from a source at a redshift Zg, would be strong evidence for 
LI breaking with 5 > 0. This is because the very large density (~400 cm~^) of 3K 
cosmic microwave photons would otherwise absorb 7-rays of energy > 100 TcV within 
a distance of ~ 10 kpc, with this critical energy reduced by a factor of ~ (1 -|- Zg)'^ for 
extragalactic sources at redshift Zg (Stecker 1969). 

3.11. Constraints on Quantum Gravity Models 

In the absence of a true and complete theory of quantum gravity, theorists have been 
suggesting and exploring models to provide experimental and observational tests of 
possible manifestations of quantum gravity phenomena. Such phenomena have usually 
been suggested to be a possible result of quantum fluctuations on the Planck scale 
Mpianck = ^Jhc/G ~ 1.12 X 10^^ GeV, corresponding to a length scale ~ 1.6 x 10~^^ m 
(Garay 1995; Amelino-Cameha, et al. 1998; Alfaro, et al. 2002). In models involving 
large extra dimensions, the energy scale at which gravity becomes strong can be much 
smaller than Mpianck, with the quantum gravity scale, Mqq, approaching the TeV scale 
(Ellis, et al. 2000, 2001). 

In many of these models Lorentz invariance can be violated at high energy. This 
results in interesting modifications of particle physics that are accesible to observational 
tests using TeV 7-ray telescopes and cosmic ray detectors. For these models, the 
constraints on S discussed in the previous section lead to significant constraints (Stecker 
2003). 

An example of such a model is a quantum gravity model with a preferred inertial 
frame given by the cosmological rest frame of the cosmic microwave background 
radiation (For an extensive discussion, see the review given by Smolin (2003). 
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In the most commonly considered of these models, the usual relativistic dispersion 
relations between energy and momentum of the photon and the electron 

E', = (8) 

E! =pI + ml (9) 

(with the "low energy" speed of light, c = 1) are modified by a leading order 
quantum space-time geometry corrections which are cubic in p ~ £^ and are suppressed 
by the quantum gravity mass scale Mqg- Following Amelino-Camelia, et al. (1998) 
and Alfaro, et al. (2002), we take the modified dispersion relations to be of the form 

- - it f^") 

E!.,l,,.l-£ (U) 

We assume that the cubic terms are the same for the photon and electron as in eqs. 
(9) and (10). More general formulations have been considered by Jacobson, Liberati 
and Mattingly (2002) and Konopka and Major (2002).+ 

As opposed to the Coleman- Glashow formalism, which involves mass dimension 
four operators in the Lagrangian and preserves power-counting renormalizablility, the 
cubic term which modifies the dispersion relations may be considered in the context 
of an effective "low energy" field theory, valid for E <^ Mqq, in which case the cubic 
term is a small perturbation involving dimension five operators whose construction is 
discussed by Meyers and Pospelov (2003). With this caveat, we can generalize the LI 
violation parameter S to an energy dependent form 

dpe dp, - Mqg 2EI Mqg' ^ ^ 

which is a valid approximation for the high energy regime E^ ^ rrif.. Note that the 
maximum velocities of particles of type i are reduced by 0{Ei/ Mqc;). 

For pair production then, with the positron and electron energy E^. c^i E,/2, 

and the threshold condition given by eq.(6) reduces to the constraint 
E^ 

M,a > ^. (14) 

Since pair production occurs for energies of at least 20 TeV, as indicated by analyses 
of the Mkn 501 and Mkn 421 spectra (De Jager and Stecker 2002; Konopelko, et al. 

+ We note that there are variants of quantum gravity and large extra dimension models which do 

not violate Lorcntz invariancc. There arc also variants for which there are only extra terms in the 
dispersion relations which are suppressed by two orders of magnitude in the quantum gravity mass and 
are therefore much smaller and do not violate the constriants given here. 



(2003), we then find the constraint on the quantum gravity scale Mqg > 0-5Mpianck 
(Stecker 2003). This constraint contradicts the predictions of some proposed quantum 
gravity models involving large extra dimensions and smaller effective Planck masses. 
Constraints on Mqg for the cubic model, obtained from limits on the energy dependent 
velocity dispersion of 7-rays for a TeV flare in Mkn 421 (Biller et al. 1999) and from 
7-ray bursts (Schaefer 1999; EUis, et al. 2003) are in the much less restrictive range 
Mqg > (5 - 7) X 10-^Mpianck. 

Within the context of a more general cubic modification of the dispersion relations 
given by eqs. (9) and (10), Jacobson, et al. (2003) have obtained an indirect limit 
on Mqg from the apparent cutoff in the synchrotron component of the in the Crab 
Nebula 7-ray emission at ~ 0.1 GeV (see Figure 10). By making reasonable assumptions 
to modify the standard synchrotron radiation formula to allow for Lorentz invariance 
violation, they have concluded that the maximium synchrotron photon energy will be 
given by -E^^max = 0.34 [eB / me){m(,/ Mqg)'"^^^ ■ This reasoning leads to the constraint 

Mqg>2.3 X l^fMpianck- 

Future observations of the Crab Nebula with the GLAST sateUite (Gamma-Ray 
Large Area Space Telescope), scheduled to be launched in 2005, will provide a better 
determination of its unpulsed 7-ray spectrum in the energy range above 30 MeV where 
the transition from the synchrotron emission component to the Compton emission 
component occurs. This will provide a more precise determination of the maximum 
electron energy in the Nebula and therefore provide a more precise constraint on the 
parameter Mqq as we have defined it here. However, this constraint will still be orders 
of magnitude above the Planck scale. 

4. Ultrahigh Energy Cosmic Ray Neutrinos 

Astronomy at the highest energies observed must be performed by studying neutrinos 
rather than photons because the universe is opaque to photons originating at rcdshift z at 
energies above 100(1 + z)'"^ TeV owing to interactions with the 2.7 K radiation (Stecker 
1969) and even lower energies owing to interactions with other sources of extragalactic 
radiation (see section 2.4.) On the other hand, the cross section for neutrino detection 
rises with energy (see, e.g., Gandhi et al. 1998) making the detection of neutrinos easier 
at higher energies. Although no ultrahigh energy cosmic neutrinos have yet been seen, 
there are various potential production mechanisms and sources which may produce these 
particles and their study can yield important new information for physics, cosmology, 
and astrophysics. 

J)..!. Neutrinos from Interactions of Ultrahigh Energy Cosmic Rays with the 3 K 
Cosmic Background Radiation 

Measurements from the COBE (Cosmic Background Explorer) convincingly proved that 
the universe is filled with radiation having the character of a near-perfect 2.725 K black 




Figure 18. The flux from photomeson production viapj2.7K followed by tt decay. 
Curve (1) is the calclated flux without redshift evolution obtained by Stecker (f979). 
The fluxes obtained by Engel et al. (2001) with redshift evolution of the proton sources 
oc (1 + z)™ with ni=3 and 4 respectively are given by curves (2) and (3). 

body, which is a remnant of the big-bang. As discussed in section 3.2, protons having 
energies above 100 EeV will interact with photons of this radiation, producing pions 
(Greisen 1966; Zatsepin and Kuz'min 1966). Ultrahigh energy neutrinos will result 
from decay of these pions (Stecker 1973, 1979). The spectrum calculated by Stecker 
(1979) without ultrahigh energy source evolution is shown in Figure 18 along with two 
flux spectra assuming source evolution oc (1 + z)"^ calculated by Engel et al. (2001). 

By extrapolating the present measurements of the flux of such high energy protons 
(see section 3), it can be shown that measurable numbers of high energy neutrinos 
can be detected using imaging optics aboard satellites looking down at the luminous 
tracks produced in the atmosphere by showers of charged particles produced when these 
neutrinos hit the nuclei of atoms in the atmosphere (see section 5). 



4-2. Neutrinos from Active Galactic Nuclei 

Quasars and other active galactic nuclei (AGN) are most powerful continuous emitters of 
energy in the known universe. These remarkable objects are fueled by the gravitational 
energy released by matter falling into a supermassive black hole at the center of the 
quasar core. The infalling matter accumulates in an accretion disk which heats up 
to temperatures high enough to emit large amounts of UV and soft X-radiation. The 
mechanism responsible for the efficient conversion of gravitational energy to observed 
luminous energy in not yet completely understood. If this conversion occurs partly 
through the acceleration of particles to relativistic energies, perhaps by the shock formed 
at the inner edge of the accretion disk (Kazanas and EUison 1986), then the interactions 
of the resulting high energy cosmic rays with the intense photon fields produced by the 
disk at the quasar cores can lead to the copious production of mesons. The subsequent 
decay of these mesons will then produce large fluxes of high energy neutrinos (Stecker 
et al. 1991; Stecker and Salamon 1996b). Since the 7-rays and high energy cosmic rays 
deep in the intense radiation field of the AGN core will lose their energy rapidly and not 
leave the source region, these AGN core sources will only be observable as high energy 
neutrino sources. 

Radio loud quasars contain jets of plasma streaming out from the vicinity of the 
black hole, in many cases with relativistic velocities approaching the speed of light. In 
a subcategory of quasars, known as blazars, these jets are pointed almost directly at us 
with their observed radiation, from radio to 7-ray wavelengths, beamed toward us (See 
sect. 2.3.) It has been found that most of these blazars actually emit the bulk of their 
energy in the high energy 7-ray range. If, as has been suggested, the 7-radiation from 
these objects is the result of interactions of relativistic nuclei (Mannheim and Biermann 
1989; Mannheim 1993), then high energy neutrinos will be produced with energy fluxes 
comparable to the 7-ray fluxes from these objects. On the other hand, if the blazar 
7-radiation is produced by purely electromagnetic processes involving only high energy 
electrons, then no neutrino flux will result. 

4.3. Neutrinos from Gamma-Ray Bursts (GRBs) 

GRBs are the most energetic transient phenomenon known in the universe. In a very 
short time of ~ 0.1 to 100 seconds, these bursts can release an energy in 7-rays alone 
of the order of 10^^ erg. They are detected at a rate of about a thousand per year 
by present instruments. It has been proposed that particles can be accelerated in these 
bursts to energies in excess of 10^° eV by relativistic shocks (Waxman 1995; Victri 1995). 

It is now known that most bursts are at cosmological distances corresponding to 
moderate redshifts (2; ~ 1). If cosmic-rays are accelerated in them to ultrahigh energies, 
interactions with 7-rays in the sources leading to the production of pions has been 
suggested as a mechanism for producing very high energy neutrinos as well (Waxman and 
Bahcall 1997; Meszaros and Waxman 2001). These neutrinos would also arrive at the 
Earth in a burst coincident with the 7-ray photons. This is particularly signiflcant since 



the ultrahigh energy cosmic rays from moderate redshifts are attenuated by interactions 
with the 2.7 K microwave radiation from the big-bang and are not expected to reach 
the Earth themselves in significant numbers (Stecker 2000). 

4-4- Neutrinos from Topological Defects and Dark Matter 

Topological defects are expected to produce very heavy particles that decay to produce 
ultrahigh energy neutrinos. The annihilation and decay of these structures is predicted 
to produce large numbers of neutrinos with energies approaching the predicted energy 
of grand unification {e.g., Bhattacharjee and Sigl 2000; see section 3.5.). The discovery 
of such a large flux of neutrinos with a hard spectrum and with energies approaching the 
energy scale predicted for grand unification would be prima facie evidence for a unified 
gauge theory of strong and electroweak interactions. 

An early inflationary phase in the history of the universe can also lead to the 
production of ultrahigh energy neutrinos. Non-thermal production of very heavy 
particles ( "wimpzillas" ) may take place. (See section 3.5.3.) These heavy particles may 
survive to the present as a part of dark matter. Their decays will produce ultrahigh 
energy particles and photons via fragmentation, including ultrahigh energy neutrinos 
(see, e.g., Barbot et al. 2003.) 

4-5. Z-hurst Neutrinos 

The observed thermal 2.7 K cosmic background radiation which permeates the universe 
as a relic of the big-bang is accompanied by a 1.96 K cosmic neutrino background of the 
same thermal big-bang origin (see e.g., Kolb and Turner 1990.) It has been proposed 
that high energy neutrinos interacting within the GZK attenuation distance with the 
copious 2 K blackbody neutrinos and annihilating at the Z-boson resonance energy can 
produce the observed "trans-GZK" air-shower events (Weiler 1982, see section 3.5.2.) 
The resulting Z-boson decays to produce a shower of leptons, photons and hadrons, 
"Z-burst". Ultrahigh energy neutrinos are produced by the decay of the tt s 
in the Z-burst. Approximately 50 neutrinos are produced per burst event; 2/3 of them 
are i/^'s. Because the annihilation process is resonant, the Z-burst energy is unique. It 
is Ez-M = 4[m^(eV)]-i ZeV. 

The Z-burst hypothesis is based on the assumption that there exists a significant 
flux of neutrinos at ~ 10 ZeV, perhaps from topological defects. Some predictive 
consequences of this hypothesis are: (a) that the direction of the air showers should be 
close to the directions of the sources, (b) that there may be multiple events coming from 
the directions of the strongest sources, and (c) that there exists a relationship between 
the maximum shower energy attainable and the terrestrially-measured neutrino mass. 

As was discussed in section 3.5.2, this production mechanism is quite speculative 
at best. 
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4-6. Neutrino Oscillations and Tau-Neutrino Observatons 

Recent observations of the disappearance of atmospheric i^^u's relative to by the 
Kamiokande group and also the zenith angle distribution of this effect (Fukuda et al. 
1998, 1999), may be interpreted as evidence of the oscillation of this weakly interacting 
neutrino state ("flavor") into another neutrino flavor, either z/^-'s or sterile neutrinos. A 
corollary of such a conclusion is that at least one neutrino state has a finite mass. This 
has very important consequences for our basic theoretical understanding of the nature 
of neutrinos and may be evidence for the grand unification of electromagnetic, weak and 
strong interactions. It is also evidence for physics beyond the "standard model" . 

If i/ju's oscillate into i/^-'s with the parameterization implied by the Super- 
Kamiokande measurements and the solar neutrino observations, (Ganzalez-Garcia and 
Nir 2002) then the fluxes of these two neutrino flavors observed from astrophysical 
sources should be equal. This is because cosmic neutrinos arrive from such large 
distances that many oscillations are expected to occur during their journey, equalizing 
the fluxes in both flavor states. Otherwise, the fluxes of z/^'s from such sources would 
be much less than those of i/^'s because u^'s are produced abundantly in the decay of 
pions which are easily produced in cosmic sources, whereas i/^-'s are not. 

Upward- moving atmospheric showers induced by ~ 100 TeV and traced back to the 
direction of a cosmic source such as an AGN or a GRB at a distance of 1 Gpc would 
occur even if the difference of the squares of the oscillating mass states were small as 
~ 10~^^ eV^. Thus, the search for upward moving showers from cosmic z/^'s, which can 
propagate thorugh the Earth through regeneration at energies above 10^"^ eV (Halzen 
and Saltzberg 1998; see section 5.2), provides a test for cosmic high energy z/^'s from 
neutrino oscillations. 

Another important signature of ultrahigh energy z/^'s is the "double bang" which 
they would produce. The first shower is produced by the original interaction which 
creates a r lepton and a hadronic shower. This is followed by the decay of the r which 
produces the second shower bang (Learned and Pakvasa 1995). The two bangs are 
separated by a distance of ~ 91.4 //m times the Lorentz factor of the r. 

4-7. Neutrino Flux Predictions 

Figure 19 illustrates some high energy neutrino fiux predictions from various 
astrophysical sources discussed above as a function of neutrino energy. Note that the 
curves show the differential + z^ flux multiplied by E^,. In this figure, z/^ <-> z/,- 
oscillations are assumed to reduce the cosmic high energy neutrino fluxes (including 
those shown in Figure 18) by a factor of 2. 

In the energy range of 0.1 to 100 PeV, the AGN neutrino flux may dominate over 
other sources. However, neutrinos from individual 7-ray bursts may be observable and 
distinguishable via their directionality and short, intense time characteristics. The time- 
averaged background flux from all bursts is shown in the figure. 

In the energy range E > 1 EeV, neutrinos are produced from photomeson 
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Figure 19. Some examples of neutrino flux predictions with <-> oscillations 
taken into account: Atmospheric and AGN fluxes (Stecker and Salamon 1996b); 
Photomeson production via pj2.7K (as in Figure 18); Topological defects (Sigl, et al. 
1998; mjf = 10^^ GeV supersymmetric model); Z Bursts (Yoshida, Sigl, and Lee 1998, 
m^ = 1 eV, primary flux oc E~^)\ 7-ray bursts (Waxman and Bahcall 1997). 

interactions of ultrahigh enrgy cosmic rays with the 2.7 K background photons. The 
highest energy neutrinos {E > 100 EeV) are presumed to arise from more the speculative 
physics of topological defects and Z-bursts. 

4-8. Observational Signatures and Rates 

4-8.1. Signatures The proposed high energy neutrino sources also have different 
signatures in terms of other observables which include coincidences with other 
observations (GRB's), anisotropy (Z-burst's), and specific relations to the number of 
hadronic or photonic air showers also induced by the phenomena (topological defects, 
Z-burst's, and 2.7K photomeson neutrinos). 

The distiguishing characteristics of astrophysical neutrino sources are summarized 
in the Table 1. In the table, the characteristic neutrino energy for photomeson processes 
is defined as ~ 1.8 x 10~^Mp/e where e is the mean energy of the photon field (Stecker 
1979), except for the case of GRB neutrinos where the energy is boosted by a factor of 
where F ~ 300 is taken to be the bulk Lorentz factor of the GRB fireball (Waxman 
and Bahcall 1997). 

The distribution of the atmospheric depth of neutrino interactions is approximately 
uniform due to the extremely long interaction path of neutrinos in the atmosphere. This 
offers a unique signature of neutrino-induced airshowers as a significant portion of the 
neutrino interactions will be deep in the atmosphere, i.e. near horizontal, and well 
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Table 1. Distinguishing characteristics of the different sources of ultra- high energy 
neutrinos (revised from Cline and Stecker 2000.) 

separated from airshowers induced by hadrons and photons. 

At ultrahigh energies, the cross sections for v and V interactions with quarks become 
equal (Gandhi et al. 1998). The interactions produce an ultrahigh energy lepton which 
carries off ~ 80% of the incident neutrino energy. The remaining 20% is in the form 
of a hadronic cascade. Charged current neutrino interactions will, on average, yield 
an UHE charged lepton and a hadronic airshower. At these energies, electrons will 
generate electronic airshowers while /x's and r's will produce airshowers with reduced 
particle densities and thus, reduced fluorescence signals. As discussed above, he z/^'- 
induced showers have a "double-bang" signature. For example, a 10 EeV, r decays after 
traveling jct — 500 km, producing a second airshower which is very well separated from 
the first, hadronic airshower at the neutrino interaction point. 

5. Present and Future Detectors 

Of the ground-based ultrahigh energy arrays, the AG AS A array of particle detectors in 
Japan is continuing to obtain data on ultrahigh energy cosmic ray-induced air showers. 
Its aperture is 200 km^sr. The HiRes array is an extension of the Fly's Eye which 
pioneered the technique of measuring the atmospheric fluorescence light in the near UV 
(300 - 400 nm range). This light is isotropically emitted by nitrogen molecules that 
are excited by the charged shower secondaries at the rate of ~4 photons per meter per 
particle. The estimated aperture of the HiRes monocular detector is ~1000 km^sr at 
100 EeV after inclusion of a 10% duty cycle (Sokolsky 1998). 

The southern hemisphere Auger array is expected to be on line in the near future. 
This will be a hybrid array which will consist of 1600 particle detector elements similar to 
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Figure 20. Two OWL satellites in low-Earth orbit observing the flourescent track of 
a giant air shower. The shaded cones illustrate the field-of-view for each satellite. 

those at Haverah Park and three or four fluorescence detectors (Zas 2001). Its expected 
aperture will be 7000 km^sr for the ground array above 10 EeV and ~ 10% of this 
number for the hybrid array. 

The next big step will be to orbit a system of space-based detectors which will 
look down on the Earth's atmosphere to detect the trails of nitrogen fluorescence light 
made by giant extensive air showers. The Orbiting Wide-angle Light collectors (OWL) 
mission is being proposed to study such showers from satellite-based platforms in low 
Earth orbit (600 - 1200 km). OWL would observe extended air showers from space via 
the air fluorescence technique, thus determining the composition, energy, and arrival 
angle distributions of the primary particles in order to deduce their origin. Operating 
from space with a wide field-of-view instrument dramatically increases the observed 
target volume, and consequently the detected air shower event rate, in comparison to 
ground based experiments. The OWL baseline configuration will yield event rates that 
are more than two orders of magnitude larger than currently operating ground-based 
experiments. The estimated aperture for a two-satellite system is 2.5 x 10^ km^sr above 
a few tens of EeV after assuming a 10% duty cycle. Figure 20 shows how two OWL 
satellite detectors will observe the track of an ultrahigh energy event from space. OWL 
will be capable of making accurate measurements of giant air shower events with high 
statistics. It may be able to detect (9(1000) giant air showers per year with E > 100 
EeV (assuming an extrapolation of the cosmic ray spectrum based upon the AGASA 
data). 

The European Space Agency is now studying the feasibility of placing such a light 
collecting detector on the International Space Station in order to develop the required 
technology to observe the fluorescent trails of giant extensive air showers, and to serve 



as a pathfinder mission for a later free flyer. This experiment has been dubbed EUSO 
(Extreme Universe Space Observatory; see the website http://www.euso-mission.org.) 
Owing to the orbit parameters and constraints of the International Space Station, the 
effective aperture for EUSO will not be as large as that of a free flyer mission. 

A recent compendium of papers on observing giant air showers from space may be 
found in Krizmanic, Ormes and Streitmatter (1998). 

5.1. Detection of Neutrino Events from Space 

The proposed Orbiting Wide-angle Light collectors satellite mission, OWL, would have 
the unmatched capacity to map the arrival directions of cosmic rays over the entire sky 
and thus to reveal the locations of strong nearby sources and large-scale anisotropies, 
this owing to the magnetic stiffness of charged particles of such high energy. As shown 
in Figure 20, this mission will involve at lease two free-fiyer light collecting satellites and 
will thus allow a stereo reconstruction of the shower tracks. With such a detector system, 
one can investigate energy spectra of any detected sources and also time correlations 
with high energy i/'s and 7-rays. 

Preliminary Monte Carlo simulations for an OWL space-based detector (J. 
Krismanic, private communication) have indicated that charged current electron 
neutrino interactions can be identified with a neutrino aperture of 20 km^-ster at a 
threshold energy of 30 EeV and the aperture size grows with energy oc E'O-^^a .^j^j^ ^j^g 
increase in neutrino cross section (Gandhi et al. 1998). Event rates can be obtained 
by convolving this neutrino aperture with neutrino flux predictions and integrating. 
Note that the neutrino interaction cross section is included in the definition of neutrino 
aperture. Assuming a 10% duty cycle of the experiment, The event rates from 
several possible UHE neutrino sources as shown in Figure 19 are found to be as 
follows: neutrinos from the interaction of UHE protons with the microwave background 
{vIi.ik) — 5 events per year; topological defects — 10 events per year; neutrinos from 
Z- bursts — 9 events per year. 

5.2. Upward Cerenkov Events from Cosmic Tau-Neutrinos 

The ensemble of charged particles in an airshower will produce a large photon signal 
from Cerenkov radiation which is strongly peaked in the forward direction and which is 
much stronger than the signal due to air fluorescence at a given energy. This translates 
into a much reduced energy threshold for observing airshowers via Cerenkov radiation. 
As this signal is highly directional, an orbitting instrument will only observe those events 
where the airshower is moving towards the experiment with the instrument located in 
the field of the narrow, Cerenkov cone. Thus, it is possible ot observe upward moving 
events from z/^'s which have propagated through the Earth (see section 4.5.) 

Virtually all particles, including neutrinos with £■ > 40 TeV, are attenuated by the 
Earth. However, i/^'s will regenerate themselves, albeit at a lower energy, due to the fact 
that both charged and neutral current interactions will have a v^- in the eventual, final 
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state (Halzen and Saltzbcrg 1998). The Vt interactions produce a leading high-energy 
r which then decays to produce another v^. 

For interactions in the Earth's crust, the r's will have a flight path of length 
7CT (~ 50 m at 1 PeV before they decay Those events which occur at a depth less than 
7cr will produce a r coming out of the Earth and generating an airshower. For a target 
area of 10^ km^, this yields a target mass of 10^(£',^(GeV) metric tons, e.g., 10^^ metric 
tons at an energy of 1 PeV for the Earth as an effective detector mass. 

Preliminary investigations of the response of an OWL space-based detector have 
indicated that the experiment would have a threshold energy of about 0.1 PeV to 
upward, Cherenkov airshowers. Assuming that the Super-Kamiokande atmospheric 
neutrino results are due to ^ Vj. oscillations, the predicted AGN flux (Stecker and 
Salamon 1996b) indicates that OWL could observe several hundred Vr events per year. 
Thus OWL would be able measure the flux of putative AGN neutrinos and observe their 
oscillations into i^^-'s. 

A downward looking balloon borne experiment to detect microwave electromagnetic 
pulses from particle cascades produced in the Antarctic ice sheet by upward moving 
ultrahigh energy neutrinos is (as of this writing) scheduled to fly in December of 
2003. It is called ANITA (Antarctic Impulsive Transient Antenna) (See the website 
http:// www. ps. uci .ed u /%7a nita / presentations, htm I . ) 
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